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In some interesting work of James Lidsey, the dynamics of Friedmann-Lemaitre- Robertson- Walker 
(FLRW) cosmology with positive curvature and a perfect fluid matter source is shown to be modeled 
in terms of a time-dependent, harmonically trapped Bose-Einstein condensate. In the present work, 
we extend this dynamic correspondence to both FLRW and Bianchi I cosmologies in arbitrary 
dimension, especially when a cosmological constant is present. 
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i 1| I. INTRODUCTION 

I . The general feature of this paper is a connection between a non-gravitational system and a gravitational system. 

r^ Such connections of course are of growing interest and importance. More specifically, extending methods initiated 

in pjj], we set up a correspondence between Bose-Einstein condensates governed by a time-dependent, harmonic 
trapping potential and both Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology and Bianchi I cosmology, in 

i__ij an arbitrary dimension, and with a non-zero cosmological constant in both cases. The correspondence is presented 
by way of two tables (Table I and Table II) that match cosmological parameters (scale factors, scalar pressure and 
energy density, Hubble parameter) with wavepacket parameters given in terms of the harmonic trapping frequency 
ui(t) and moments Ijit),j = 2,3,4 (with hit) > 0) of the wavefunction of the Gross-Pitacvskii equation - under the 

1^-. ' assumption that the atomic interaction parameter is a constant in time. Here t is "laboratory time" that one passes 

r»o , to from cosmic time r in the Einstein field equations. The moments satisfy the conservation law 

^, 

|^! . X(t) = 21 2 {t)h(t) - if(t)/4 = a constant A. (1) 

o, 

Moreover, l2(t) gives rise to a solution X(t) =' I 2 (£) of the classical Ermakov-Milnc-Pinncy (EMP) equation [1^ 



- w+ ^it)X=^. (2) 



On the other hand, one knows also that the Einstein equations for a FLRW universe, and even for the anisotropic 
Bianchi I and Bianchi V universes, admit a formulation in terms of a suitable EMP equation, of classical or of 
generalized type [51481 fill Il5l [l6l|. This information coupled with equation (2) allows one to proceed in setting up the 
desired correspondence - at least in the FLRW and Bianchi I cases. 

For certain cosmological models, for example a stiff fluid model, the frequency u(t) can be determined, and thus the 
external potential is explicated. A new issue arises however when a non-vanishing cosmological constant is present. 
Our analysis shows that in this case one must employ elliptic functions to solve the appropriate moment equations 
that arise - a matter discussed in section V. 

As indicated in [12J, through the condensed matter - cosmology correspondence via EMP equations (as considered 
here) there is the increased possibility for insight into the hidden symmetries of these systems. Further connections 
between non-gravitational and gravitational systems are explored in [J], for example. 

II. AN EMP FORMULATION OF FLRW AND BIANCHI I 

We begin by formulating the d— dimensional FLRW and Bianchi I equations with a cosmological constant A^ (for 
d > 3) as a single EMP equation of classical type; we do this (for the record) as it has not been done before in 



arbitrary dimensions, even though the extension is straightforward. The classical EMP equation is of the form 



Y + A(t)Y = 



JL 

Y 3 



(3) 



for Y = ^r, and for some constant /i; compare equation (2). 

We consider a scalar field <p and a potential V, and let a(r) denote the scale factor, in which case the FLRW 

equations assume the following form, for the Hubble parameter H = ' — : 



H l + — 



2A„ 



2K d 



a 2 (d-l)(d-2) (d-l)(d-2) 



A'\2 



r'Y D 



(4) 



<t>'4>" + {d- i)H{4>') 2 + (V o 4>)4>' = o 



(5) 



with k = 0, — 1, or 1 the curvature parameter, Kd — SuGa for Gd = the gravitational constant and A^ the cosmo- 
logical constant. 

Suppose /(£) > is a function with inverse function T(t) (i.e. f(T(r)) = r, T(f(t)) = i) such that T'(t) = a(r). 
Define 



K(t) = / -a(/(t)),0i(t)^ / X/(t)). 

The method of [1JJ shows immediately that equations (4) and (5) lead to the EMP equation 



Y(t) 



nK d 



2(d-2) 



^(i) 2 F(i) 



F(t)a 



(6) 



(7) 



For the anisotropic d— dimensional Bianchi I cosmological model with metric ds 2 = —dr 2 + X\(r) 2 dx\ + 
Xd-iir) 2 dx 2 d _ 11 the field equations take the form 



^HiH k =K d 



Kk 



+ V0(j) 



+ A d , 



(8) 



J2 {Hi + Hf) + J2 H l H * = -K d 



if-' 



Kk 

l,kfi 



- Vo 



A, 



(9) 



where i,l,k £ {1, . . . , d— 1} and Hi = ' Xi/Xi. By making the substitution Xi(t) = i?(r)e Q ^ T ^ for functions R(t) > 
and <xi{t) satisfying ai(r) + • • • + cvd_i(r) = 0, one sees that the field equations (8) and (9) can be written as 



(d-l)(d-2) 



Hi- 



DKd 

R2(d-i) 



= K r 



r 



+ V0(j} 



A d , 



(10) 



(d-2)H' R 



(d-l)(d-2) 



Hi 



DKd 

R 2{d-1) 



-K A 



- Vo, 



Ad 



(11) 



since R = {X 1 X 2 ■ ■ ■ X d -i) 1/{d - 1 \ where H r (t) d =' R'(t)/R(t), D d = S=mr d J2i<k( H i ~ H kf ■ D is a constant 
quantity by a simple lemma which states that for any differentiable function g(r), any positive diffcrcntiable function 
R(t), and M <E M, the function g(r)R(r) M is constant if and only if g / (r) + Mff(r) „/V = 0. By equating the left-hand 
sides of any two Einstein equations (9), indexed by i ^ j, we see that the latter equation holds for g = Hi — Hj and 
M = (d - 1), so that (Hi - H^R^ 1 is constant for all i,je{l,...,d- 1}. 

Similar to the above argument, by taking f(t) > to be the inverse of T(t) such that T'(t) = R^)^^ 1 ^ , and by 

defining Y(t) d =' R(f (t)f d ~ l \ ^(t) d = 4>(f(t)), we can show that the Bianchi I field equations (10), (11) lead to the 
classical EMP 



VM j. Mj^ i w a vw -2(d-l)K d D 
Y(t) + -Jd~2T Mt) Y{t) = (d-2)Y(tf ■ 



(12) 



III. ELLIPTIC FUNCTIONS INTERLUDE 



In the case of a non-vanishing cosmological constant, that we shall give attention to, the following differential 
equation (which is of some independent interest) 



4 y2 



(13) 



arises, where A > 0,C ^ 0; sec equation (35). Its solutions involve elliptic functions. Before setting up the dynamic 
correspondence between cosmological and condensate systems we present a brief interlude regarding equation (13), 
which in particular will allow, conveniently, for the introduction of some notation needed later. 

We shall need the elliptic functions EF(x,k), E(u,k) of the first and second kind, respectively, with modulus k 
II I MM given by 



EF{x,k) de d- f 

JQ 



(It 



E(u,k) d = 



o v/(l-i 2 )(l-fc 2 i 2 ) 

sn(u,k) 

dn 2 (v, k)dv 



(14) 



1 - fc 2 t 2 

1-t 2 



dt. 



Sometimes, as is usual, the modulus k is suppressed in the notation and one writes E{u), sn(u), for example, for 
E(u,k), sn(u,k). Given that C / 0, we can define constants m = m(C),n = n(B,C) and Weierstrass invariants 
g 2 =g 2 (B,C),g 3 =g 3 (A,B t C),p = p(B,C),q = q(A,B,C) by 



def. C i lo „ 4 



// 



def. 



\c\ 

D 



l ■ 3 def. 4 . 

( i-e- rn 6 = 77 ), 



.93 = 

def. 



3C" 

def. 2B 3 



27C 2 

93 . 
4 



!J2 

-2A 



def. mB 2 



3C* 



mnB, 



(15) 



B :i 



54C 2 



def. 32 

p = ~T = 

A 

f 2- 



mB 2 
' 12(7' 



The point of these definitions is that if u(£) is the function given by v(t) 
to the differential equation 

( TOt; + ^)V ^ 4 
4(?j 3 + pw + g) to 2 ' 



de/. 



[y(t) — n]/m, then equation (13) transforms 



(16) 



!>■{■ 



as a direct computation reveals. Moreover, if w(x) is an inverse function of y(t) then z{x) =' w(mx + n) is an inverse 
function of v(t) so that from equation (16) the equation 



z'{xf 



m 2 (mx + ^) 



16 X(x) 
is directly derived. This means that we can focus on solving 

mB 1 



,X(x) d ^x 3 



px 



z\x) = ± 



which involves the elliptic integrals 



dx 



l2C^XJx) 
xdx 



(17) 



(18) 



4 ^X(x)_ 
that in turn involve a consideration of the roots of the 



^X(X) J y/X(xj^ 

cubic equation X(x) = for their evaluation. The case of interest here is when X(x) = has one real root T\ and 
two complex roots r2,r 3 - the other cases being simpler to deal with. Necessarily r2 and r 3 are complex conjugates: 
7"3 = T2- The condition that X{x) = indeed has a single real root r\ and two complex conjugate roots r2,r 3 is that 

its discriminant A = -4p 3 - 27q 2 d = i [g| - 27 g 2 } d = ±^4 [2B 3 - 27 AC 2 ] (since m 3 = £) should be negative: 
2B 3 — 27 AC 2 < (given that A > 0), which we therefore assume. 



Associated with the roots are useful parameters a, p, Q,t\,ti 



def. T def. ,~ def. 

a = Imr 2 , p = -r x /2, g = 



vw 



def. 



ri + V '(p - r i) 2 + ° 2 = n + VV + ct 2 , 



(19) 



ti = ' r\ - \f{p-ri) 2 +<r 2 = n - V9/0 2 



Here <t 7^ 0, since otherwise r2 would be a second real root. The cubic X{x) admits the factorization X{x) = 

(a* — ri)[(x — p) 2 + cr 2 }, which shows that (since ir^O) X(x) > for x > r\. Thus x > r\ will be a convenient 

assumption in our discussion of functions like X(x) and z'(x) in (18), for example, where , 1 appears. The elliptic 

modulus fc that will be employed in (14) is given by 



k"'d- 



\ 



3ri 

2 



= +\ 



<)rr 



y/9p 2 + g 2 + 3p 
2V9p 2 + cr 2 



(20) 



We note that < fc < 1. For ct ^ 0, ^/VTct 2 > ^/V = 3|p| > ±3p =■> 2^9p 2 + cr 2 = ^9p 2 + a 2 + ^9p 2 + a 2 > 
v/9/9 2 + o- 2 + 3p > so that 1 > 



^/9p 2 +<T 2 +3p 



> 0: i.e. l>fc 2 >0=^l>fc>0, as desired. 



With the preceding notation and definitions in place, we can now construct a crucial function u(x) = u(x, fc) that 
facilitates the expression of solutions of equation (18). For 



e{x) de J- i 



{x-hf 



(x-t 2 ) 2 

[2-c-(*i+t 2 )](*i-*2) 
(x - 1 2 ) 2 



n < x, 



and for EF(x, k) in (14) and k in (20), we set 



def 



(x) = u(x, k) 2- EF{y%x), k). 



(21) 



(22) 



Here t 2 < 7*i < x => x - t 2 ^ in (21). Also 6(ti) = 1, but for r\ < x ^ t\ (i.e. (x - ti) 2 > 0) 9(x) < 1. Moreover, 
2x > 2ri = h + 1 2 (by (19)) =*> 0(x) > (since ti > t 2 ). That is, n < x ^ t\ =» < 9(x) < 1 =*> < V^M < 1 in 
(22). Since sn(y) is the inverse function of EF(x), we can provide a second description of u(x). Namely, y/0(x) = 

sn(EF(y/6{x))) = sn(u(x)), and therefore cn(u(x)) = \J\ — sn 2 (u(x)) = \J\ — 9{x) = ' J^I^I (since (again) x > t 2 
for x > r*i ) =>■ 



t(x, fc) = en' 



\x-h\ 

x - t 2 



r\ < x. 



(23) 



In terms of the elliptic functions EF(x,k),u(x,k) defined in (14) and (22) (or (23)), again with the elliptic modulus 
fc specified in definition (20), equation (18) is solved as follows, for an integration constant zq: 



z{x) = ± 



m 



2 gh , mgB\ 



12C 



J u ( x ) 



m 2 , , .. m 2 \/X(x) 
—EF(u(x)) + V )> 
2.g 2{x - t 2 ) 



zq 



(24) 



for j*i < x < ti, and 



z(x) = ± 



m m gt 2 mgB 



! v^R 



12Cj u { x ) + -EF(u(x)), ^^^ 



zo 



(25) 



for t\ < x. m,g,ti,t 2 arc defined in (15), (19), and we assume that A > in (13), and 27 AC 2 > 2B 3 so that n is 
the unique real root of the cubic equation X(x) = 0, for X(x) in (17) with p, q there also defined in (15). 



def. def. 

Going back to the definition z(x) = w(mx + n) (for n = B/3C in (15)), where w(x) is an inverse function of 
y{t), we sec that the initial differential equation (13) is solved implicitly by way of w(x) given by 

w{x) = z {^r) ' (26) 

for z(x) given by (24), or by (25). 

IV. THE CORRESPONDENCES 

The Lidsey correspondence between BEC's and cosmology originates by way of a comparison of equation (2) with 
equation (7), resulting in the following table. 

TABLE I: BEC o FLRW correspondence 



h <* a 2 

h ^ 2(aH) 

ll/4J 2 o H 2 

[(d-l)(d-2)h-A d ]/K d o P<t> 

■2)w 2 h~(d-l){d-2)h + A d ]/K d o P4> 



Similarly, the correspondence in the case of the Bianchi I cosmology originates by comparing equations (2) and (12) 
for which one obtains the following table. 

TABLE II: BEC o Bianchi I correspondence 



h ++ 2(d-l)(R^ d -^H R ) 

Ii/ih ♦+ (d - \) 2 H 2 R 

[{fc§/4 - A„] /K d o P4> 

[^^ 2 h-^h + A d ]/K d «. pt 

A <-> -2(d-l)K d D/(d-2) 



SOME EXAMPLES 



The moments are known to satisfy the equations [12t [1 

h(t) = 0, I 2 (t)=h(t), (27) 

I 3 (t) = -2u(t) 2 I 2 (t)+4I 4 (t), j 4 (t) = =^-J 3 (t). 

At this point we assume an equation of state p^ = (7 — l)p^>, with 7 > 0. We indicate how to extend the discussion 
in section IV of [12J]- The initial step is to equate the fifth BEC entry in Table I with (7 — 1) times the forth entry 
there and solve for ui 2 . The result is that 



which by (27) gives (since ^3 = I2) 



h = 



7(d-l)/ 4 7 A d 



2 h 2(d - 2)/ 2 



(29) 



(30) 



(29) being a first order, linear differential equation consequently has the solution 

j = _£ A d 

4 j7(d-l)/2 "•" (d_l)(d-2)' 

for an integration constant a, which plugged into equation (28) yields the Ad— independent result 

W - .[ 7 (d-l)+2]/2- ^1) 

J 2 

By equations (1) and (27) 

A = 2I 2 h - l'i/4. = 2I 2 h - ll/4, (32) 

which with the help of equation (30) can be re-written as 



U 2a 2A d I 2 



I [~ f (d-l)-2]/2 (d-l)(d-2) 



- A. (33) 



Note that the above assumption that 7 > (so in particular, 7 7^ 0) rules out the un-wanted conclusion u) = 0, 
by equation (31) (or (28)), and also the conclusion that Ii(t) is a constant function, by equation (29). Equations 
(31) and (33) govern the time-dependent trapping frequency and thus the external potential V(r,t) = ui(t) 2 r 2 /2 also, 
although equation (33) is a bit complicated. If A^ = 0, for example, it has the implicit solution 



9/ 2+1 T? flj-1 1.3 1 1 , A/g 

2 2-ri \ 2 ^ q ^ 2' 2 ^ q' 2a . 

±t + t Q (34) 



2a(q + 2) 

for q = [7(0? — 1) — 2]/2, and for an integration constant to- 

Consider the choice 7 = 6/(d — 1), for example, which corresponds to a stiff perfect fluid. Then equation (33) 
becomes 



1$ = 2a_ 2A d 7 2 

4 1$ + (d-l)(d-2) 



^ = ^-A+ „ "7 oV (35) 



This is equation (13) for A d =' a, B d =' A, and C d =' 2A<j/(d— l)(d— 2), with A d 7^ 0, where we note that ^4 = a > by 
(30), since 7 > 0. Equation (35) can therefore be solved (implicitly) in terms of the elliptic functions EF(x, k), E(u, k) 
in definition (14). Namely, the inverse function of ^(i) is given by w(x) = z (^p-) where z{x) is given by equation 
(24) or (25), according to equation (26) and the notation of section III. The hypothesis there is that 27 AC 2 > 2B 3 . 
This means that 54oA 2 ! /(d — l) 2 (rf — 2) 2 > A 3 , which as we have seen corresponds to the condition of a negative 
discriminant A. The cases A > and A = are treated in the Appendix. By equation (31) 

w = 74 • (36) 

1 2 



In case Ad = 0, one can write equation (33) as \ I2.&I2. 1 \J la. — XI 2 = ±di, which integrated gives —-^^/2a — \I\ = 
±t + t , or 

I 2 2 {t) = ^--A\[±t + t a f, (37) 

for an integration constant t . By equation (37), equation (36) is explicated. 

Another choice of interest is 7 = -At , (for a universe dominated by matter when d = 4) in which case equation 

(33) reads If/4 = A/I 2 -B + CI 2 , or 



dl 2 = ±2t + 1 (38) 



y/CI$ -Bh+A 



where A = 2a, B = A, C — 2Ad/(d — l)(d — 2), and to = an integration constant. For a general value of C the integral 
in (38) can be expressed explicitly in terms of the elliptic functions EF(x,k),E(u,k) of the first and second kind in 



definition (14) by use of Maple, for example. Thus again for A^ 7^ 0, I 2 (t) can be determined implicitly, and moreover, 
by equation (31), 

^ = £ (39) 

In the particular (easier) case when A^ = (i.e. C = 0), for example, the integral in (38) is an elementary function. 
Namely, equation (38) reduces to the equation 



J A- BI 2 Jh A / / BI 2 \ , , 

- B + B^ arCtan [\Ia^bT 2 ) = ±2t + ^ (40) 

again for A =' 2a > 0, B = A, say for < I^t) < -§ = x> wnere already iaCO > an d ( as we have seen) a > 0. 



Thus one needs that A > 0; in [12| the choice A = 1 is made. 

As a final example, regarding table I, we take 7 = ^rj- Then equation (33) can be written as 



lTdh 



2a + CI 3 / 2 - XI 



1/2 



±2t + t 0l (41) 



def 

again for C =' 2Ad/(d — l)(d — 2). However, the integral here is non-tractable unless C = (i.e. A^ = 0), in which 
case its evaluation gives the implicit equation 



6a 2 / y/XI, 

—rj^arctan \ - y 2a — XL 

ka - XI 2 /2 



U 2 1/2 



3a/ 2 1/4 , IT 
A 2 A 



= ±2i + t (42) 



for I 2 (t). Also for 7 = 3/(d - 1), to 2 = ia/ll' 2 by equation (31). 

To close things out we present a few examples regarding Table II for a Bianchi I cosmology, where we maintain the 
equation of state p<f, — (7 — l)p0,7 > 0. In place of equations (28) and (29), one quickly checks that the equations 
lu 2 = 7I4//2 - l{d - l)A d /{d - 2)I 2 ,h + {l/1)hh/h = (7/2) (d - l)A d i 2 /(d - 2)I 2 follow by Table II. The latter 
equation has solution 

a ( d — 1\ . , , 



'2 



by which the former equation can be written as 



J 2 

which again is A^-independent, and where (again) a is an integration constant. By equations (32) and (43) we deduce 
that 



is the Bianchi I version of equation (33). The equation 

dl 2 



2a/ 2 1_7/2 + bI 2 -X 



±2t + 1 (46) 



def 

is a re-expression of equation (45) for b =' 2(d — VjAd/yd — 2). 

An obvious solution for I 2 is obtained by choosing 7 = 2, for example, in which case equation (46) reads (for 
A d ± 0) (2/b)y / bI 2 (t)-X + 2a = ±2t + t a , or 



Equation (44) then assumes the explicit form u! 2 (t) = 2ajl\(t). 

As a second example, choose 7 = 1, A^ = (i.e. 6 = 0). Then equations (44) and (45) read 



3^ 



±2t + t 



, 2 (t)=a// 3/2 (i), 



(48) 



since 



/(it I I fi f 

— '" = A\ a\fx — c (a\fx+2c)l'Sa 2 . More generally, the integral / —========= fo 

xfaJx - c v J v / ax 1 /(»+i) - c 



_.rn = 0,l,2,3--- , 

y/ay/x- c v J VaxVl'i+i) _ c 

which normally involves the hypergeometric function 2-P1 , can be explicitly computed. Thus for A^ = 0, we obtain a 
family of examples by choosing 7 = y n ='2(1— 7^ l,n = 0,l,2,3.... For the record, we arc able to derive the 
general formula for the corresponding integral in (46), 



d.r 



Vax 1 /(«+i) _ " a n + 



2(n + 1) / — 777— 7; ^ n!(ax 1 /(™+ 1 ) - C )"-V 

— -\/ax 1 /( n + 1 ) — r > — 

1 ^i!(n-i)!(2n-2j + l) 



(49) 



i=o 



The case n = 1, was just treated. If n = 2, for example, 7 = 73 = 4/3 and the corresponding integral in equation (46) 
is computed by the formula 



dx 



2Vax 1 /3 _ c 



3a 2 a; 2/3 + 4acx 1/3 + 8c 5 



vWVS-T^ 5a 3 

If n = 3, then 7 = 73 = 3/2 and the integral in equation (46) is computed by 

VaxV* - c [5a 3 x 3 / 4 + 6a 2 x 1 / 2 c + Sao^V + 16c 3 



y/axV* - c 35a 4 
In general for 7 = 7„ , equation (44) shows that 



2na 



(n + l)f 2 - 1/(n+1) 



(50) 



(51) 



(52) 



APPENDIX: COMPUTATION OF SOME ELLIPTIC INTEGRALS 

The problem of solving equation (13) has been reduced to that of solving equation (18), which in turn is a matter 
of computation of the elliptic integrals Ij(x) = J 7== ■ For the reader's convenience we provide the result, which 



y/X(w) 



can be deduced from formulas in [3j coupled with a few extra arguments. We use freely the notation of section III. 

There are three cases: A < 0, A > 0, and A = 0. First assume that A < 0. That is, 27 AC 2 > 2B 3 so that X(x) = 
has a single real root r\. Then, omitting integration constants, we have 



h{x) 



gu(x) for n < x < t\ 
—gu(x) for ti < x 



(Al) 



h(x) 



(«(*) - E(u(x)) + ^) 

(u(x)-E(u(x))-^ 



+ t%gu{x) for ri < x <ti 
— t2gu{x) for ii < x 



for u(x) — u(x, k) in definition (22). 



If A > 0, i.e. 27 AC 2 < 2B 2 , then X(x) = has throe distinct real roots a, b, c, say a > b > c. In this case we now 
define k, u(x) by 



k d g- , /*Z£ u{x) = u(a!j fc) */■ BJ r (\/*-^l, fc I . (A2) 



a — c \\ x — b 

Then X{x) = (x — a)(x — b)(x — c) and for x > a 

2u(x) 



Io(x) 



\fa— c 



(A3) 

I\(x) = 2\J a — c [dn u(x)tn u(x) — E(u(x))] H , 

Va- c 

where tn a; = ' snx/cnx as usual. 

If A = (the final case), then X(x) = has at least two real roots: X(x) = (x — a) 2 {x — c) for real numbers a, c. 
Thus the Ij(x) are elementary functions computable by a calculus table of integrals, depending on whether c = a or 
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